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q . Abstract 

o 

kjjy Event-by-Event fluctuations for high energy heavy-ion collisions are investi- 

gated. The ^-measures are calculated for pions. It is found that multiparticle 
\ symmetrization correlation has great influence on the value of if the phase 

Cf~) ' space density is high. Previous derivation of <£-measure has neglected multi- 

particle quantum correlation. So the explanation of ^-measure as a measure 
■ of dynamical correlation is reasonable only when multiparticles quantum cor- 

relation effects are smaller. We argue that the ^-measure might be applicable 
■^-j- ! at the SPS energy where the pion phase space density is much less than one. 

o 
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One of the main goals of relativistic heavy-ion physics is the discovery of the new state 
of strongly interacting matter that is called the quark-gluon plasma(QGP) []J. As one can 



[ not detect QGP directly, lots of indirect signals are suggested It has been pointed out 
that the study of event-by-event fluctuations in high energy heavy-ion collisions may provide 
us more information about the experiment • For examples, event-by-event fluctuations 
may provide us information about the heat capacity [|,|ff|], possible equilibration of the 
system |9|-|rT|j or about the phase transition ||T8|,|8|j. Experimental group NA49 has already 



studied the ^-measure in Pb+Pb collisions at SPS energy |19| . Recently Bialas, Koch |2( 
and Belkacem et. al. |HJ have found a relationship between the moments of event-by-event 
fluctuations and the inclusive correlations. In this paper we will extend the derivation given 



in [2I|, and connect the result of Ref. [20| with Ref. [21]. It has been suggested in Ref. || 
that if there is no special dynamics involved in nuclei and nuclei (AA) collisions, then the <£>- 
measure for AA should be the same as the ^-measure for proton and proton (pp) collisions. 
So by comparing ^-measures for pp and AA, one can get information on the dynamics that 
is involved in AA collisions. However, the derivation in Ref. || has neglected the multi- 
particles Bose-Einstein symmetrization effects among the identical particles. The difference 
between ^-measure for AA and that for pp collisions may be due to the multi-particles Bose- 
Einstein correlation so that the statement of Ref. 11 given above is no longer valid anymore. 



In Ref. [12| , the author has assumed that the pions produced from AA collisions has a Bose- 
Einstein distribution form in momentum space and has used this spectrum distribution to 
calculate ^-measure and found that ^-measure of Bose-Einstein symmetrization is around 



the prediction of the experiment. But it has been shown in Ref. 22] that only when the 



1 



volume goes to infinity, the spectrum distribution of pions can be of the function form of 
Bose-Einstein. In this paper, we will calculate the ^-measure of pions in heavy-ion collisions 
by using a pion distribution function in both coordinate and momentum space. As has been 



shown in Ref. |23j that when coordinate space volume goes to infinity, multi-particle Bose- 
Einstein correlation will make the pion spectrum distribution change to Bose-Einstein form 
in momentum space. So the calculation in Ref. |L2] can be achieved in our calculation 
by letting the volume goes to infinity. On the other side, the effects of the finite size of 
AA system on ^-measure can be studied in our calculation. This paper is arranged in the 
following way: In section II, the relationship between ^-measure and particles inclusive 
distributions are derived. In section III, we first derive the formulas which will be used for 
the calculation of ^-measure of pions, then we calculate the ^-measure for different system 
size and temperature. For comparison, the ^-measure of pions based on the method given 



in Ref. |L2[ are also calculated. In section IV, a simple model is given which demonstrates 
explicitly that the derivation given in Ref. [17] has neglected multi-particles Bose-Einstein 
symmetrization effects and conclusions given in Ref. [|17j only holds when the multi-particles 
symmetrization effects are smaller. Finally, Conclusions are given in Section V. 



II. THE RELATIONSHIP BETWEEN ^-MEASURE AND PARTICLES 

INCLUSIVE DISTRIBUTION 



The quantity Z(y) for an arbitrary event is defined as P.|20 



N 



Z(y) = 5>(pi) 



(1) 



i=l 



here N is the multiplicity of the event and y(pi) is the function of momentum pi in the 
event [] 

Then the k-th order moments of Z can be written down as [pO|,pl[ 



x If y(pi) is total momentum, y(pi) = \pt\. If y(pi) is transverse momentum, y(pi) = \Jp\ + p|> If 
y(pi) is the rapidity, then we need to know the energy of particles i. If y(pi) is electric charge of 
the particles, we also need to know the electrical charge of particles i. So Eq.(|l|) can be generalized 
to 

N 

Z(y)=5>(Pi,Si), (2) 
i=i 

Here Si represents other property of particles i, like mass and charge. One can simply write Eq.(|2|) 
as 

N 

Z(y) = J2vi (3) 
i=i 

Here is the function of i-th particle. 
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M 



(4) 



i=l 



where M is the total number of events and Nj is the multiplicity of the j-th event. 

In the following we will use iV-body distribution function f(N, pi,- • -,Pn) which is 
normalized as 



/ dpi • • • dp N f(N, Pl , • • •, PN ) = P(N). 
Here P{N) is the normalized multiplicity distribution. 



(5) 



Following Ref. ||21|| , we define the reduced s-body distribution functions for a system of 
indistinguishable N— particles as: 



f (Af=7)T / dp s+ i • • dp N f(N,p ir ■ -,p N ) s < N 
f 8 (N, pi, Ps ) = I N\f(N, pi, • • •, p N ) s = N 

I s > N. 

But we should point out that this definition is different from that in Ref. [f2l|[ p| 
According to Eq.(0), One can write out the fc-order moments of Z as |2"0|.|2~T|| : 



(6) 



-, k 



(Z(y) k }=Y.[dpi---d P K Ez/( P i) /(iV,pi,...,p N ). 

N=0 J Li=l 

To be clear, the first to fifth order moments of Z can be expressed as: 

oo - 

(z(y)) = Y, d Pi 2/(Pi) A(^Pi)> 

N=0 

00 r r r 

(Z(y) 2 ) = E J d Pi ?/ 2 (Pi) /i(^Pi) + J dpidpi y(pi)y(p a ) / 2 (iV,pi,p 2 ) 



7V=0 



(9) 



(10) 



OO r „ „ 

(^%) 3 > = E y d Pi y 3 (pi) A(iv,pi) + 3 y d P idp 2 y 2 (piMp 2 ) / 2 (iv,pi,p 2 ) 



V=0 



+ y dp 1 dp 2 dp 3 y(pi)y(pa)3/(Pa) /s(^Pi,P2,Pa) 



(12) 



2 In Ref. [21], the authors have defined a function f s (N, pi, ...,p s ) as 

AT! 



f s (N,p 1: ...,p s 



dp s+1 ...dp^f N (N, pi, ... P n)- 



(7) 



(JV-s)! 

It is clear that when s = N, the above formula is "bad" defined while this situation does not exist 
in our definition (see Eq.@). On the other hand, one can derive from Eq.(^) of this paper that 



f s (N, pi, p s ) = - y dp a+ x...dp N f N (N, pi, ...pn)- 
From Eq.(fJ) and Eq.(||), we know that this two definitions (Ref. |2l| and here) are different. 



(8) 
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(Z(y) 4 ) = £ f /dpi 2/ 4 (pi) /i(JV sPl ) + 4 / d Pl d P2 y 3 {p 1 )y(p 2 ) / 2 (iV,p 1)P2 ) 
+3 / dp 1 dp 2 Z/ 2 (pi)y 2 (P2)/2(A^,pi,p 2 ) 

+ 6 / rfpirfp 2 rfp 3 y 2 (Pl)Z/(P2)y(P3)/3(^, Pl,P2,P3) 

+ / rfpirfp 2 rfp 3 rfp 4 y(pi)y(P2)y(P3)y(p4) /4(-W,pi,P2,P3,P4) 



(Z(yf) = E [/ dpi y 5 (pi) /i(JV, Pi) + 5 / rf Pl rfp 2 y 4 (piM P 2) / 2 (iV, pi, Pa ) 
+ 10 1 rf Pl rfp 2 y 3 (pi)y 2 (p 2 )/ 2 (iV, Pl)P2 ) 

+ 10 /"rfpirfp 2 rfp 3 l/ 3 (pi)y(p 2 )y(p3)/3(iV,Pl,p 2 ,P3) 

+ 15/" dpidp 2 dp 3 2/ 2 (Pl)y 2 (p 2 )y(P3)/3(^,Pl,p 2 ,p 3 ) 

+ 10/" dpirfp 2 rfp 3 ?/ 2 (pi)y(p 2 )y(p3)y(p4)/4(iV,pi,p 2 ,p3,p4) 

rfpirfpadpa^rfps 2/(PiMP2MP3MP4MPs) fs (JV, Pi , p 2 , P3 , P4, Ps) 



The s particles inclusive distribution can be written down as 



^=0 



Then we have p0 



(13) 



• (14) 



p(Pi, P2, -, Ps) = E fs( N > Pi> Pa> ■••> Ps) (15) 

7V=0 

with 

- 00 AT' 

y dpi...dp.p( Pl , Ps ) = e p ( N ) (N ^ s)l = ( N ( N - - s + !)>■ ( ie ) 



(z(y)> = / y(p)p(p)dp, (17) 

(^(y) 2 ) = y l/ 2 (p)p( P )d P + y dpidp a j/( Pl )y(pa)p(pi,Pa), (18) 

<^) 3 > = y Z/ 3 (p)p( P )d P 

+3 y rf Pl d P2 y 2 ( Pl )y( P2 )p( Pl , P2 ) 

+ y dpidpatip32/(pi)3/(pa)3/(ps)/?(pi,P2,Ps) J (19) 
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(Z(y) A ) = J y 4 (p)p(p)dp + 4 J dp 1 dp 2 y 3 (p 1 )y(p 2 ) 

+3 J dp 1 dp 2 y 2 (p 1 )y 2 (p 2 )p(p 1 ,p 2 ) 

+6 J dp 1 dp 2 dp s y 2 (p 1 )y(p 2 )y(p 3 )p(p 1 ,p 2 ,p 3 ) 

+ / dpidp 2 dp3dp42/(pi)y(p2)y(p3)y(p4)p(pi,P2,P3,P4) 



(20) 



and so on. One can use Eq.flD or Eqs(0,0,|3) to calculate (Z(y) K ) separately to examine 
the consistency between those equations. The discussions about the consistency between 
those equations can be found in Ref. |2D . 



III. ^-MEASURE OF PIONS WITH MULTI-PARTICLES BOSE-EINSTEIN 
SYMMETRIZATION CORRELATIONS 

In the following we will only consider pions(7r + or tt~), the abundant particles produced 
in high energy collisions. The N pions distribution function can be written down as p3| 



with 



/(iV, Pl ,..,p N ) = P(iV) 



Pv(Pl, -,Pn) = ^P (iV) (Pl ) P<x(l))--P (iV) (PN,P ( T(N)) 



Pv(Pl,P2, -,Pn) 
J rfpi...rfp N P7v(pi, p N ) 



(21) 



(22) 



Here p^(pi, Pj) is the Fourier transformation of the pion source distribution function 
^^(x, K), and g( N \x., K) is the probability of finding a pion at point x with momen- 
tum K. The superscript iV denotes those quantities are calculated for a system with pion 
multiplicity N. a(i) denotes the i-th element of a permutations of the sequence {1, 2, • • -, s}, 
and the sum over a denotes the sum over all s\ permutations of this sequence. 
Using Eqs. (|,|1],|2D, we have HUf: 



f s (N, pi, - • p 



i=s mi=l m,2=l 



wW(i\T) 

i— mi— rrt2 m s _2— 1 

) P CT (l))G m V 2 ) (P2,P ff (2) 

m s _i=l cr 

■ • ■ G {N) 



■ in s -x IPs— 1) Pcr(s-l))G i _ mi nis 



,( Ps)P(J ( S) VW(iv-i) 



(23) 



with 



u (N) (N) = ^J dpi...dp N iMpi,...,PN) 

= ^E- w (at-,)/^(p,p), 



j=l 

W W(0) = l,Gf = ^(p,p) 



(24) 
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and 

Gf^q) = / P W (P,Pi)W iV) (Pi ! P2) ■ ■ ■ dpi-iP W (Pi-i,q)- (25) 

The superscript iV of and G^(p, q) denote o/"^ and (p, q) are calculated using 
the quantities p^(pi,pj). 

The s-pion inclusive distribution can be expressed as 

oo 

p(Pi, ••;Pb)= E f»( N > Pi' -•' P g ) 

Ar=0 

oo AT i— (s— 1) i— mi— (s— 2) i~m 1 -rri2 m s _2-l JV / AT _ ,-\ 

N=s i=s mi=l m 2 =l m 3 _i=l V J V <r 
G r ^ ) (PbP CT (i)) G? ^ ) (P2 ) P<r(2)) ' ' ' G m a -\ (Ps— 1) P<r(s— 1) —mi m s _i (Ps,Pa( S )) 

uo^ N \N — m\ — ni2 — ... — m s ) 



E---E E m)- 

mi=l m s =l A r =mi+m2+...+m s 



»P<t(1)) ■ • ■ G £!(PnP<r(B))- (26) 

Neglecting symmetrization effects, that is p^ N \phPj) = 0, we have 

J N \N) = ± (27) 

and 

G'f ) (Pi,P j ) = 0, A;^lU Pi ^ Pj . (28) 

Thus Eq.(|23|) reads 



f s (N, Pl ,..., Ps ) = -^L G f( Pl , Pl )...Gf( Ps , Ps VW(iV- s ) 

= P (^)^rZ^)j^(Pi' Pi)-^(Ps> Ps)- (29) 

From Eq.([0|) or Eq. (|26l) we have 

00 

p( Pl , Ps ) = E N ( N -l)...(N-s + l)P(N)p N (p 1 , p 1 )...p 7V ( Ps , Ps ). (30) 

N=s 

If we neglect multi-particle symmetrization(or antisymmetrization) effects and the correla- 
tion between particle spectrum distribution and multiplicity, we immediately come to the 
conclusion that p iV (p,p) = p(p)/(N). Of course this conclusion is based on the fact that 
we have neglected the dynamical correlation and kinematic correlation. 
In calculation, one normally choose the variable y(p) as [f| 

y(p) = x(p) - W, x = J dpx(p)j^-. (31) 
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If we neglect all kinds of correlations in the system, the first three moments of Z read: 

(Z k ) = J dpy k (p)p(p) (l<k< 3). (32) 



For high-order moments, Eq.(|32|) does not hold anymore. This can be easily seen from 
Eqs. (^3| . [T^ . [20| ). This result is different from Ref. |2l| where Eq.(|32]) have been thought to 
be valid for all-order moments. Mrowczynski has also noticed this fact in Ref. |I7|| . So 
^-measure 



l/k 



(v 



(2<k<3) 



with 



y 



(N) 



p(p)y k (p)dp, 



(33) 



(34) 



is related to the quantum many-body correlations, the correlation between the multiplicity 
and momentum distribution, and the dynamical correlations, the last one has already been 
discussed in Ref. ||. 

As we have no interest in the first two kinds of correlations, we will eliminate those 
correlations from the ^-measure. One may get rid of multiplicity and momentum correlation 
by redefining the variable y(p) as 



yW(p) = s(p) 



x 



(TV) 



with 



x 



(N) 



dpx(p) 



N 



(35) 



(36) 



Here ■^ l( ^' p ' > is the normalized particle spectrum distribution of pions with multiplicity N. 
After this definition one may get rid of the possible correlation between the multiplicity and 
momentum. In the actual calculation, one may calculate the average value of x(p) for each 
event(denotes as x e ), then uses the redefined y(p) = x(p) — x e and Eq.(^) to calculate (Z k ). 
For the y k , one may calculate it using following equation 



1 M N i 

2/ fc = T7EE%(Pi) fc - 

1V1 3=1 i=l 



(37) 



Here M is the total number of events and Nj is the multiplicity of the j-th event. yj(p) = 
x(p) — Xj is the quantity defined for the j-th event, and Xj is mean value of x(p) for the 
j-th event. Using this method, one can calculate the $- measure. 

In the following we will calculate the quantum correlation effects on ^-measure. As- 
suming (/^(x, K) as a Gaussian function in both coordinate space and momentum space, 



P (Pi)Pj) can De expressed as P|P5 



P (7V) (Pi,Pj: 



9\x, 
1 



Pi + Pj 



exp(z(pi — pj) ■ x)dx 



(2ttA 2 



Texp( _ (Pl -P^W )exp( 



vPi + Pj, 

8A(Ny 



(38) 
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Here R{N) and A(iV) are the source radius and temperature for the events with multiplicity 
N. Of course the source radius and temperature may be different for different event, even 
though those events have the same multiplicity. But at this calculation we will neglect this 
fluctuation. In Fig.l, $2(Pt) vs - mean pion multiplicity are shown Q It is interesting to 
note that $ value will increase as the pion multiplicity goes up. For the same multiplicity, 
$ value will decrease as the volume or temperature drops. All these can be understood 
in the following way: $2(^4) will become bigger as the phase space density (~< N > 
I (RA) 3 ) becomes higher. Of course, the $ value should also depend on the pion multiplicity 
distribution p3| , |25| , |27|] , but this effect will be neglected at present calculation. Here we would 
like to point out that many authors have noticed that Bose-Einstein symmetrization may 
affect the value of ^-measure. For example, In Ref. [21|, similar formulas as Eqs (|lO , |Tl|, 12) 



of the paper have been derived, but the functional form of p(pi,...,p s ) have never been 
presented. On the other hand, Mrowczynski has studied Bose-Einstein correlation effects 
on the ^-measure by assuming the particle spectrum distribution is Bose-Einstein form 



Here, we will explain how to calculate the ^-measure given in Fig.l. According to Eq.(P,B3), 
we have 

* 2 (ft) = (|^) 1/2 - (^r / p[p)? t dv)V\ (39) 



(N) ' ^ (N) 

with 



Pt=Pt~ (pt), ipt) = T^y / p(p)Pi • ( 10) 



Here pt is transverse momentum (pt = ypf+p^), (N) is the mean multiplicity. According to 
Eq.fll8|), we have 



(Z 2 (p t )) = J Ptp(p)dp + J dp 1 dp 2 pitP2tp(pi,P2)- (41) 
According to Eq,(p6|), the single and two-particles inclusive distribution can be written down as 



00 00 



p(p)= E E m)^ u ^G mi (p, P ), (42) 



mi=l N=rrii 



oo 00 



\ - \ - ^ , ,u(N - mi - 1712) 

p(pi, P2 )= E E £ p w — — 

mi=l mi=2 Af=mi+m2 v ' 

[G , mi (pi,p 2 )G m2 (p 2 ,pi) + G m 1 (Pl> Pl)C m2 

(P2,P2)]- (43) 

In the calculation we assume that R and A are the same for all multiplicity, so we neglect the super- 
script N in Eq.(|26|), From the function p^ N \pi,Pj) given in Eq.([38|), one can calculate G mi (pi,p2) 
according to Eq.(p5|) and uj(N) according to Eq,(p4|). Choosing P(N) as Poisson form, that is 



P(N) = exp(— (N)) , we can calculate p(pi,P2) and p(p). Bring all those formulas in eq.(| 
we get $ 2 (pt)- 
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which corresponds to R is infinity for our calculation | 23| . However Neither Belkacem et. 
al. nor Mrowczynski has calculated multi-body Bose-Einstein symmetrization effects on 
$— measure as done in this paper. For comparison the results of Ref. [0 are also shown 
in Fig. 10. The specific volume (V = L 3 = 1000/m 3 ) and temperature (T = 150MeV) are 
chosen for the formulas of Ref. ]12j so that (r 2 = x 2 + y 2 + z 2 ) boXtL=10 f m ~ (r 2 ) 

Gaussian, R=3fm 

and (p 2 )r=i50MeV ~ (p 2 ) A =o.2GeV- It is interesting to notice that our results is larger than 
the results of Ref. |E| when the phase space (~ ((r 2 ) (p 2 )) 3 ^ 2 ) are the same. $2(^4) for Ref. 
| 12|| will also become bigger if the phase space density becomes higher. Here we should point 
out that for (N) ~ 100, which is the mean multiplicity of 7r + or 7r~ measured in SPS energy, 
our Q(pt) is much larger than the experimental value of 4.6 ± 1.5 ||19|| . However for large 
volume (R = 5.3 fm and A = 0.2GeV) or high temperature (R = 3fm and A = 0.355GeV), 
we find that $ ~ 7MeV for < >~ 100 which is consistent with the $ value calculated by 
two-particles correlation (5 ± 1.5Mel^) for Pb+Pb at SPS |19 |. Taking V = L 3 
and T = 150MeV, we find that $ value of Ref. fTM 



(16) 3 fm 3 



is around 6MeV which is consistent 
with our results and data. This suggests that for large volume (R = 5.1 fm or L = 16 fm 



which is around the HBT radius of NA49 p8|), high-order Bose-Einstein symmetrization 



effects are not significant for SPS energy. So our calculation is consistent with the results of 
Ref. H. 



4 Here we will explain how to calculate the ^-measure using the method given in Ref. JT2|. Ac- 
cording to Ref. [12 1, we can infer that 

V 1 

P(P) = (27r) 3 exp(/3(£-/j))-l' (44) 

Here V = L 3 is the volume, (3 = T~ l is the inverse temperature, fi is the chemical potential. 
E = \/ m 2 + p 2 with m being the particles mass and p its momentum. The mean multiplicity (N) 
can be expressed as 

(AO = J dpp(p). (45) 



The mean {pt) reads 



(«> = S -&$±. (46) 



It is easily checked that 

\2 exp(l3(E-u)) 



{ z 2 (p t) ) idpip t -(p t )) 2 ^ m^^ _ (47) 



Bring those formulas in Eq.(p^), we can get <&2(pt)- 
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IV. THE SHORTCOMING OF PREVIOUS DERIVATION OF ^-MEASURE 



In the following we show that the previous derivation of the ^-measure is based on a 
model which has no quantum correlation among particles. We assume that in AA collisions, 
particles are produced by each nucleon nucleon (NN) collisions and the number of collisions 
of NN has a distribution function P(M). Here M is the number of collisions in each event 
and the collision center has a distribution of p(x). The probability amplitude of finding a 
pion at point x can be written down as [^9j 

M M 

J(x)=^J ? ;(x)=5:j(x-X i ) (48) 
i=l i=l 



Here Xi is the collision center coordinate. In the last step of Eq.(f4*8|) we assume that Jj(x) 
is only the function of x — x;. The probability amplitude in the momentum space can be 



expressed as 29 



M , M 



J(p) = / c?xJ(x) exp(ip • x) = E / exp(ip ■ x) J(x — Xi)cfx = E J(p) exp(ip • x;) (49) 
^ i=i ^ i=i 

with 

J(p) = J J(x) exp(ip ■ x)dx. (50) 
Then the s-pion inclusive distribution reads 

p(Pl,-,Ps) = E P ( M ) / rfx 1 ...ciXM^(x 1 )...^(xM)|J(p 1 )| 2 ...|J(Ps)| 2 

M M 
|E ex P(«Pl •X i )| 2 ...|E eX P( i Ps- X i)| 2 ( 51 ) 



i=l i=l 



the second line of Eq.(^) comes from Eq.(^9|) which is due to quantum effect. So the 
distribution P(X) in Ref. []I7| does not hold anymore. However if we consider classical 



particles, we immediately can get an equation similar to the expression of P(X) in Ref. [17 
As the terms in the second line of Eq. ( 51 ) change to 

M M 

IE ex p( ? pj • x 0l 2 -> E I ex p(^j • x i)l 2 = M ( 52 ) 

1=1 8=1 



we have 



M=oo 

p(Pi,p 2 ,-,Ps)= E P(M)M s \J(Pi)\ 2 -\APs)\ 2 - (53) 

A/=0 



Then the derivation in Ref [17] can be applied here. In passing, we would like to note 
that the P(X) distribution in Ref. |T7| is the source distribution, while, p(pi, ...,p s ) in our 
work represents particle inclusive distribution. For illustration, we will write out the mean 
multiplicity (N) as: 
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(iV> = J p( Pl )d Pl 

= (M) J d Pl \ J( Pl )| 2 + (M(M - 1)) J d Pl \J( Pl \ 2 \g(Pi)\ 2 (54) 

with 

fif(Pi) = / dx£f(x) exp(ipx • x). (55) 

Here / <ipi|J(pi)| 2 = no is the mean number of pions that are emitted from each NN 
collisions. The previous relation like (N) = (M)no does not hold anymore as the quantum 
effect has extra contribution to the average value of (N). Of course the last term of Eq. 
([54]) depends on the distribution function of the source, g(x), and it will become smaller if 
the collision region is larger. This can be explained in quantum mechanical context: If the 
source of AA collisions is very larger and dilute pion gas exists, quantum effect will become 
smaller. 



V. DISCUSSIONS AND CONCLUSIONS 

It has been proved that ^-measure will be significantly influenced by the multiparticle 
correlations effects if the particle phase space density is very high. However one may exclude 
the quantum symmetrization or antisymmetrization effects by cuts. One may do HBT 
analysis and determine the HBT momentum region first, then throw away the identical 
particle pairs which are affected greatly by particles symmetrization effects from each event. 
After this is done one may use the remaining particles to re-calculate the $ measure. It is 
interesting to mention the following fact: In Ref. H50 , |5I |, the authors have already calculated 



the average phase space density for pions at SPS energy, and have found that the pion phase 
space density is 0.2, which is much less than 1. Therefore multipion symmetrization effect 
may not be significant at SPS energy. For RHIC energy, it is estimated that the pion phase 
space density will become higher, so one should be careful when using ^-measure in RHIC 
energy. 

In the paper, we calculate the ^-measure for a specific source distribution function. 
There is a possibility that the results of ^-measure depends on the function form of the 
source. As has been stated in the paper that ^-measure also depend on the fluctuation 
of the source radius and temperature for different multiplicity. All this kind of work has 
not been studied in detail in the paper as the main interest of present paper is to study 
the Bose-Einstein correlation effects on the ^-measure. We find that ^-measure is greatly 
affected by Bose-Einstein symmetrization if the source density of pions is high. 

In this paper, we derive the relationship between ^-measure and particles inclusive dis- 
tributions. It is shown that besides dynamical correlation, ^-measure is also related to 
the many-body quantum (anti) symmetrization correlations. Starting from a pion source 
distribution in coordinate and momentum space, we calculate the ^-measure for pions. It 
is found that ^-measure is strongly affected by multiparticle symmetrization effects if the 
phase space density is very high. So the previous statement that ^-measure reflects the 
dynamical mechanism that exists in the AA does not hold anymore if the pion phase space 
density is high. For a very large system, the density of pions may be smaller, in this case 
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the multiparticles symmetrization effect can be neglected. Finally we argue that one may 
more or less exclude the quantum symmetrization effect by cutting small momentum pairs 
of pions from each event. 
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FIG. 1. 3>2(pt) vs - mean pion multiplicity (N). In the calculation we choose A = 0.2GeV, 
R = 3/m(dashed line) and 4/m (long dashed line) respectively. For comparison, &2(pt) for 
A = 0.25GeV, R = 3/m(dot-dashed line) are also shown in Fig.l. We choose the pion multi- 
plicity distribution, P(N), as Poisson form. The results from Ref.[ll] are also shown in Fig.l. 
We choose the volume of the pion system (box) as V = L 3 = 1000/m 3 and the temperature 
T = 0.l5GeV (dotted line) and T = 0.18GeV (solid line) respectively. 
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